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Correlation Effects in Quantum Dot Wave Function Imaging 
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We demonstrate that in semiconductor quantum dots wave functions, as imaged by local tunneling spectroscopies 
like STM, show characteristic signatures of electron-electron Coulomb correlation. We predict that such images 
correspond to "quasi-particle" wave functions which cannot be computed by standard mean-field techniques in the 
strongly correlated regime. From the configuration-interaction solution of the few-particle problem for prototype 
dots, we find that quasi-particle wavefunction images may display signatures of Wigner crystallization. 
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1. Introduction 

Present single-electron tunneling spectroscopies in 
semiconductor quantum dots 1-3 (QDs) provide spectac- 
ular images of QD wave functions. The measured in- 
tensities are generally identified with the density of car- 
rier states at the resonant tunneling (Fermi) energy, re- 
solved in cither real 4-6 or reciprocal 7-9 space. However, 
Coulomb blockade phenomena and strong inter-carrier 
correlation, which are the fingerprints of QD physics, 
complicate the above simple picture. 10 Indeed, QDs can 
be strongly interacting objects with a completely discrete 
energy spectrum, which in turn depends on the number 
of electrons, 1,3 N. Therefore, orbitals can be ill-defined, 
losing their meaning due to interaction. Also, it is un- 
clear how many electrons one should take into account 
to calculate the local density of states, as a particle tun- 
nels into a QD and the number of electrons filling the 
dot Actuates between N — 1 and N (such fluctuation is 
the origin of either the Coulomb current peak or the ca- 
pacitive signal). 9, 11-13 

Here we clarify the physical quantities actually 
probed by scanning tunneling microscopies 4-6 (STM) or 
magneto-tunneling spectroscopies 7-9, 11 of QDs, and how 
they depend on interactions. If only one many-body state 
is probed at a time, then the signal is proportional to the 
probability density of the quasi-particle (QP) being in- 
jected into the interacting QD. We demonstrate that the 
QP density dramatically depends on the strength of cor- 
relation inside the dot, and it strongly deviates from the 
common mean-field (density functional theory, Hartree- 
Fock) picture in physically relevant regimes. 

2. Theory of Quasi-Particle Imaging 

The imaging experiments, in their essence, measure 
quantities directly proportional to the probability for 
transfer of an electron through a barrier, from an emit- 
ter, where electrons fill in a Fermi sea, to a dot, with 
completely discrete energy spectrum. In multi-terminal 
setups one can neglect the role of electrodes other than 
the emitter, to a first approximation. The measured 
quantity can be the current, 4,7 the differential conduc- 
tance, 5,6,8, 12 or the QD capacitance, 9,11,13 while the 
emitter can be the STM tip, 4-6 or a n-doped GaAs con- 
tact, 7-9, 11-13 and the barrier can be the vacuum 4-6 as 

*E-mail address: rontani@unimore.it 



well as a AlGaAs spacer. 7-9,11-13 

According to the seminal paper by Bardeen, the 
transition probability (at zero temperature) is given by 
the expression jK) \M\ n(ef), where M. is the matrix 
element and n{ef) is the energy density of the final QD 
states. The common wisdom would predict the probabil- 
ity to be proportional to the total density of QD states at 
the resonant tunneling energy, e j , possibly space-resolved 
since M. would depend on the resonant QD orbital. 15 Let 
us now assume that: (i) Electrons in the emitter do not 
interact and their energy levels form a continuum, (ii) 
Electrons from the emitter access through the barrier a 
single QD at a sharp resonant energy, corresponding to 
a well defined interacting QD state, (iii) The xy and z 
motions of electrons are separable, the z axis being par- 
allel to the tunneling direction, (iv) Electrons in the QD 
all occupy the same confined orbital along z, XQd( z )- 
Then one can show 10 that the matrix element A4 may 
be factorized as 



M cx TM, 



(1) 



where T is a purely single-particle matrix element while 
the integral M contains the whole correlation physics. 

The former term is proportional to the current density 
evaluated at any point z>, ar in the barrier: 



T = 



2m* 



Xe(*0 



dXQp(^) 
dz 



Xqb(z) 



gx|(f) 
dz 



(2) 

where xe{z) is the resonating emitter state along z 
evanescent in the barrier and m* is the electron effective 
mass. The term (2) contains the information regarding 
the overlap between emitter and QD orbital tails in the 
barrier, Xb( z ) an( i XQd(-z), respectively. Since T is sub- 
stantially independent from both N and xy location, its 
value is irrelevant in the present context. 

On the other hand, the in-plane matrix element M 
conveys the information related to correlation effects: 



M 



<>v{Q) <PQD(g)dg, 



(3) 



where <£qd{q) is the QP wavefunction of the interacting 
QD system: 16 

m n(g) = (N-l\if( e )\N). (4) 
Here </>e is the in-plane part of the emitter resonant or- 
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bital, ^(f?) is the fcrmionic field operator destroying an 
electron at position g = (x, y), \N — 1) and \N) are the 
QD interacting ground states with N — 1 and N elec- 
trons, respectively (see also Sec. 3). We omit spin indices 
for the sake of simplicity. 

Results (3-4) are the key for predicting wave function 
images both in real and reciprocal space. In STM, </>e(£>) 
is the localized tip wave function; if we ideally assume 
it point-like and located at Qo, 15 i-e. 4>e(q) ~ 3(q — Qa), 
then the signal intensity is proportional to |c/?qd(£>o)| 2 , 
which is the usual result of the one-electron theory, 6, 15 
provided the ill-defined QD orbital is replaced by the 
QP wavefunction unambiguously defined by Eq. (4). 
In magneto-tunneling spectroscopy, the emitter in-plane 
wavefunction is a plane wave, 0e(#) = c lk ' e , and the 
matrix element (3) is the Fourier transform of <^qd, 
M = (fQY)(k). Again, we generalize the standard one- 
electron result 8 by substituting <pQn(k) for the QD or- 
bital. Note that M is the relevant quantity also for in- 
tensities in spacc-intcgratcd spectroscopies probing the 
QD addition energy spectrum. 12 ' 13 

3. Two-dimensional Quantum Dot 

3.1 The Non- Interacting Case 

We now apply the theory of Sec. 2 to a two-dimensional 
parabolic QD with a few strongly interacting electrons. 
The harmonic potential was proven to be an excellent 
description of several experimental traps. 3 The non- 
interacting effective-mass Hamiltonian of the i-th elec- 
tron is 

Ho(i) = ^ + ( 5 ) 

The eigenstates <fi a (o) of (5) are known as Fock-Darwin 
orbitals. 17 Their peculiar shell structure, with constant 
energy spacing hwo, is represented in Fig. 1 up to the 
third shell. 



N = 5 



N=6 



d 



Fig. 1. Electronic configuration of the non-interacting quantum 
dot ground state as the electron number, N, fluctuates between 
5 and 6. The arrows represent electrons (with their spin) filling 
Fock-Darwin orbitals. The letters identify different energy shells 
(the third-shell central orbital has s character). 



What is (Pqb{q) in the non- interacting case? Let us 
consider e.g. A = 6. Then, the 5-electron ground state 
| A — 1) appearing in the definition (4) is naturally ob- 
tained from the Aufbau principle of atomic physics: 12 
all and only the lowest-energy Fock-Darwin orbitals are 



filled in with electrons according to Pauli exclusion prin- 
ciple (left panel of Fig. 1). The 6-electron ground state 
|6) = Cp||5) is obviously obtained from |5) by adding 
a spin-up electron into the lowest-energy p-type empty 
orbital (right panel of Fig. 1); here is the perti- 
nent creation operator. By expanding the field oper- 
ator ^ on the Fock-Darwin orbital basis, &(g, s z ) — 
Y.aa fa(Q) £<r(s z ) K<x(s z ) is the spin part of the elec- 
tron wave function with eigenvalue a =f, j and spin co- 
ordinate s z ], we derive that the orbital part of the QP 
wave function is simply <£>qd((?) = <Pp{o)- 

The above result is a sensible one: as we inject e.g. via 
the STM tip an additional electron to the 5-electron 
ground state and N oscillates between 5 and 6 (Fig. I), 
the non-interacting wave function of the tunneling elec- 
tron can be regarded alternatively either as the lowest- 
energy unoccupied orbital when N = 5 (Fig. 1 left panel) 
or as the highest-energy occupied orbital when N = 6 
(Fig. 1 right panel). In the section below we consider the 
effects of electron-electron interaction. 

3. 2 Configuration- Interaction Approach to the Interact- 
ing Problem 

The fully interacting Hamiltonian is the sum of single- 
particle terms (5) plus the Coulomb term: 



N 1 

tf = 5>o(o + 5 E 



i=i 



2i£K\Qi-Qj\ 



(6) 



where k is the static relative dielectric constant of the 
host semiconductor. We solve numerically the few-body 
problem of Eq. (6), for the ground state at different 
iV's, by means of the configuration interaction (CI) 
method, 17-19 where \N) is expanded in a series of Slater 
determinants built by filling in a truncated set of Fock- 
Darwin orbitals with N electrons, and consistently with 
symmetry constraints. From the solution of the result- 
ing large-size matrix-diagonalization problem, we obtain 
eigenvalues and eigenvectors of ground- and first excited- 
states. 20 Then, we evaluate the matrix element (4), by 
decomposing \N) and — 1) on the Slater determinant 
basis: the resulting <£qd(£?) is now a mixture of differ- 
ent Fock-Darwin orbitals, with weights controlled by the 
strength of correlation. 

3.3 Tuning the Strength of Correlation 

A way of artificially tuning the strength of Coulomb 
correlation in QDs is to dilute the electron density. While 
the kinetic energy term scales as r~ 2 , r s being the pa- 
rameter measuring the average distance between elec- 
trons, the Coulomb energy scales as rj 1 . Therefore, at 
low enough density, electrons pass from a "liquid" phase, 
where low-energy motion is equally controlled by kinetic 
and Coulomb energy, to a "crystallized" phase, remines- 
cent of the Wigner crystal in the bulk, where electrons 
are localized in space and arrange themselves in a geo- 
metrically ordered configuration such that electrostatic 
repulsion is minimized. 3 In the latter regime Coulomb 
correlation severely mixes many different Slater determi- 
nants, and the CI approach is the ideal tool to quantita- 
tively predict correlation effects. 17 
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The typical QD lateral extension is given by the char- 
acteristic dot radius £qb = (h/m*^) 1 ^ 2 , £qd being 
the mean square root of g on the Fock-Darwin lowest- 
energy level ip s . As we keep N fixed and increase ^qd, 
the Coulomb-to- kinetic energy ratio A = £qd / [ajg = 
ft 2 k I '(m* e 2 ) is the effective Bohr radius of the dot] 21 in- 
creases as well, driving the system into the "Wigner" 
regime. 22 As a rough indication, consider that for Aw2 
or lower the electronic ground state is liquid, while above 
A w 4 electrons form a "crystallized" phase. 21 

4. From high to low electron density: Wigner 
crystallization 

Figure 2 shows the square modulus of the QP wave 
function, corresponding to the injection of the 6-th elec- 
tron, in the xy plane for three different values of A. As 
A increases (from top to bottom), the density decreases 
going from the non-interacting limit (Fig. 2, top panel, 
A = 0.5), deep into the Wigner regime (Fig. 2, bottom 
panel, A = 10). At high density (A — 0.5, approximately 
corresponding to the electron density n e — 3.8 x 10 12 
cm~ 2 ) the wave function substantially coincides with the 
non-interacting Fock-Darwin p orbital f p (g) of Fig. 1. By 
increasing the QD radius (and A) , the QP wave function 
weight moves towards larger values of g. By measuring 
lengths in units of ^qd, as it is done in Fig. 2, this trivial 
effect should be totally compensated. However, we see in 
the middle panel of Fig. 2 (A = 4, n e s» 1.5 X 10 10 cm~ 2 ) 
that the now much stronger correlation is responsible for 
an unexpected weight reorganization, which is related to 
the formation of an outer "ring" of crystallized electrons 
in the Wigner molecule. 21 Such tendency is clearly con- 
firmed at even lower densities (A = 10, n e w 1.3 x 10 9 
cm' 2 ). Now, together with the outer ring, a new struc- 
ture is visible close to the QD center (if g — ► then 
ifQP due to the orbital p symmetry). Such complex 
shape is consistent with the onset of a solid phase with 5 
electrons sitting at the apices of a regular pentagon plus 
one electron at the center. 21,23,24 

In Fig. 2 the absolute QP weight has been arbitrarily 
renormalized, which we believe to be a sensible procedure 
in view of comparison with experimental images. In order 
to illustrate a second correlation effect, in addition to 
shape changes, in Fig. 3 we plot the absolute value of the 
QP wave function as a function of x at y = 0, for the same 
values of A as in Fig. 2. 10 Figure 3 clearly demonstrates 
a dramatic weight loss as A is increased: the stronger the 
correlation, the more effective the orthogonality between 
interacting states. Note also in Fig. 3 that the shoulder 
of the outer QP peak close to the QD center is clearly 
visible for A = 10. 

5. Present Status of Imaging Experiments 

Among the existing imaging tunneling experi- 
ments, 4 ~ 6,8 ' 9, 11 two have specifically focused on results 
in the presence of several carriers in the QD. 9, 11 

A first experiment 9 concerned electrons in InAs self- 
assembled QDs in the non interacting high-density limit 
(A w 0.5). This work demonstrates the experimental 
imaging of an Aufbau-like filling sequence as up to six 
electrons are sequentially injected into the QDs. In par- 




Fig. 2. Square modulus of the quasi-particle wave function in the 
quantum dot plane for three different values of the dimension- 
less parameter A. This quantity is proportional to the STM signal 
when the electron number N fluctuates between 5 and 6. As A in- 
creases (from top to bottom) the density decreases and the wave 
function shape evolves from a characteristic p-type Fock-Darwin 
orbital (top panel, A = 0.5) into a complex figure peculiar of the 
"crystallized" phase (bottom panel, A = 10). The wave function 
normalization is arbitrary and the length unit is the character- 
istic dot radius ^qd- 

ticular, the specific filling sequence can be understood 
in terms of the consecutive filling of the s orbital first, 
then one of the twop orbitals of the second shell (Fig. 1), 
and eventually the other one. The above sequence differs 
from that expected according to Hund's rule, namely the 
third and fourth electrons should separately fill in the two 
p orbitals with parallel spins. 12 The deviation from the 
above rule is attributed to either piezoelectric effects or 
to a slight elongation of the InAs island shape. 9 However, 
many QDs were probed at once and the above interpre- 
tation cannot be regarded as definitive. 

In a second experiment the imaging of InAs QD hole 
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Fig. 3. Quasi-particle wave function vs. x (y = 0) for differ- 
ent A-values. We use the same parameters as in Fig. 2 except 
that now the wave function normalization is absolute. As A in- 
creases the total weight decreases from 0.97 (A = 0.5) up to 0.15 
(A = 10). 



wave functions was addressed. 11 The authors observe an 
anomalous filling sequence up to 6 holes (s, s, p, p, d, d) 
and interpret it in terms of a generalized Hund's rule for 
the two p and the two d orbitals together, namely the to- 
tal spin should be maximized as N increases as an effect 
of strong Coulomb correlation. Assuming reasonable hole 
parameters as k = 12.4, Hlvq = 25 meV, m* = 0.3m e , we 
estimate, within the simple parabolic potential model, 
the key parameter A to be 1.46. Such value is compa- 
rable to those of typical devices 12 showing "standard" 
Aufbau physics and it seems to us too small to support 
claims of qualitatively new correlation effects. 17 An alter- 
native explanation of the filling sequence could be related 
to merely single particle effects arising from the complex 
hole band structure. Specifically, assuming orbital ener- 
gies to be ordered as s, p, d and no shell degeneracy, then 
electrons should consecutively fill in such orbitals and 
Hund's rule would never hold. Such interpretation seems 
to be confirmed by independent theoretical work. 25 Fur- 
ther measurements, as a function of the magnetic field 
parallel to z, could be useful to clarify the question. 

From the above discussion it appears that experimen- 
tal investigations of QP wave functions in regimes where 
correlation effects are significant are lacking so far. We 
hope that our results will stimulate further experiments. 

6. Conclusions 

In conclusion, wc have shown that QP wave functions 
of QDs are extremely sensitive to electron-electron corre- 
lation, and may differ from single-particle states in phys- 
ically relevant cases. This result is of interest to predict 
the real- and reciprocal-space wave function images ob- 
tained by tunneling spectroscopies, as well as the inten- 
sities of addition spectra of QDs. We believe that our 
findings will be important also for other strongly con- 
fined systems, like e.g. nanostructures at surfaces. 26 
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